Quantum effects on one-dimensional collision dynamics of fermion clusters 
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Recently, many experiments with cold atomic gases have been conducted from interest in the non- 
equilibrium dynamics of correlated quantum systems. Of these experiments, the mixing dynamics 
of fermion clusters motivates us to research cluster-cluster collision dynamics in one-dimensional 
Fermi systems. We adopt the one-dimensional Fermi-Hubbard model and apply the time-dependent 
density matrix renormalization group method. We simulate collisions between two fermion clusters 
of spin-up and spin-down, and calculate reflectance of the clusters R changing the particle number 
in each cluster and the interaction strength between two fermions with up and down spins. We 
also evaluate the quasi-classical (independent collision) reflectance R'^'^ to compare it with R. The 
quasi-classical picture is quantitatively valid in the limit of weak interaction, but it is not valid when 
interaction is strong. 



INTRODUCTION 

In recent years non-equilibrium dynamics of cold atom 
systems has emerged as interesting problems, because 
the cold atom systems are ideally isolated systems vifhich 
can be directly explored experimentally. Recent rapid 
progress in manipulation techniques has enabled experi- 
mentalists to modify strength and sign of interaction be- 
tween cold atoms by Feshbach resonance [1], to make a 
spin-dependent trapping potential to trap the cold atoms 
of each spin separately, and to change the trapping po- 
tential suddenly to explore the dynamics of quantum 
quench. 

Of these experiments, the dynamics of the mixing of 
two spin components of fermion gas [2] motivates us to 
study collision dynamics between two fermionic atom 
clusters in one-dimensional Fermi systems. The fermion 
cluster-cluster collision dynamics is so simple to compare 
it with classical dynamics, but the quantum effects such 
as Fermi statistics on collision dynamics have not been 
clarified. Furthermore, cluster collision dynamics is one 
of the basic concepts of dynamics, which can be applied 
to diverse quantum dynamics. Therefore the research of 
fermion cluster collision dynamics brings a new perspec- 
tive of quantum dynamics. 
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FIG. 1. (a) Trapping potentials of the two types of fermions 
and particle density of the trapped fermion clusters {n — 6), 
and (b) the shared potential. 
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THEORETICAL CONDITIONS 

Simrdations are conducted in the one-dimensional sys- 
tem which consists of two-component fermions. In the 
system, at first, the fermions of each component are 
trapped by separate harmonic potentials. The trapping 
harmonic potentials have the same shape (the cycle of 
motion of a fermion in the potential is T), but they are 
spatially separated (FigjIJa)). The particle density has 
some peaks (the number of peaks is the same as the par- 
ticle number) and has Gaussian-like tails. 

Then, we suddenly change the trapping potentials into 
a new shared potential (we set this time i = 0). The 
new shared potential has the same shape as before, but 



its bottom is midway between the bottoms of the two 
initial trapping potentials (FiglTJb)). When the inter- 
action between the two fermions is zero, each fermion 
cluster moves to the opposite side in the half cycle T/2, 
without changing the shape of the particle number den- 
sity. On the other hand, if the interaction is finite, some 
particles are reflected and return to the initial side, and 
the others pass through and move to the opposite side 
(FiglHa)(b)). At i = T/2, the particle number density is 
nearly zero again at the center of the potential, so we can 
obtain the number of particles which return to the initial 
side; this number is denoted by N^^'^^ [n is the number 
of particles in the initial cluster). Therefore we can ob- 
tain the reflectance Rn = N^^^ /n, and the transmittance 
r„ = 1 — Rn ■ Clearly, if the interaction between the two 



fermions is zero, i?„ = 0, and if the interaction is infinite, 
Rn = 1 because the system is one-dimensional. 
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FIG. 2. (a) Particle density at f = r/2 and the interaction 
u = 9.05[ft77r"^] (n = 6), and (b) particle density at i = r/2 
and interaction u — 144.8[?!.77r~^] (n — &). 



SIMULATIONS BY T-DMRG METHOD 

We adopt the one-dimensional Fermi-Hubbard model 
and apply the time-dependent density matrix renormal- 
ization group(t-DMRG) method to simulate the system. 
We take TVgitc = 4i:) + l sites numbered 0, 1, ..., 2L», ...,41? 
(D is an integer), and the system size is L = 4a; the site 
is at a; = —2a, the site 2D is at a; = 0, and the site 41? is 
at a; = 2a. The lattice constant is Aa; = ajD. The trap- 
ping potential is V'^*^{x) = ^muj^{x + a)^ for the spin- 
down particles, and V^^*^{x) 



mu) [x—a) for the spin- 
up particles, where m is the particle mass, and lo = I-k jT 
and T is the cycle of the oscillation. The continuum limit 
is recovered in the limit of A'sitc ^ oo. The edge of the 
particle density decays in about r] = 1/ ^/moj. At i = 0, 
we change the trapping potentials to V^^{x) = ir/iw^a;^. 
The Hamiltonian is 

i i,a 

Here J = {2m Ax^)^^ is the hopping constant, and U is 
the on-site interaction between a spin-up fermion and a 
spin-down fermion. In the continuum limit, the interac- 
tion term reproduces the contact interaction u5{xi —X2), 
in which u = UAx is the interaction strength. 



Here we show the results for the system size L = 20ri 
(a = Sry) and D = 50. The system size is large enough 
compared to the clusters' expanse, so the particles do 
not escape from the system, and we assume that D is 
large enough so that the effect of the discretization is 
insignificant. 



RESULTS 
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Interaction u [tin/T] 

FIG. 3. (a) Reflectance Rn for u = UAx between 0.200 and 
2.260 [hriT''^], and (b) ratio of reflectance Rn/Ri- 



Figure ^a) shows the reflectance i?„ for the interac- 
tion u between 0.200 and 2.260 [hriT-'^]. The simple 
theoretical calculation predicts R\ (x v?. The figure im- 
plies Rn oc V? for any n. We plot the ratio of i?„ to i?i in 
FiglHJb), and it is observed that i?„ ~ nRi in the limit 
of u — ; 
u^O. 



0. Therefore the results indicate i?„ oc nu as 



On the other hand, Fig|4l^a) shows the transmittance 
T„ for u between 144.8 and 3276.8 [firfT-^]. Theoreti- 
cally, Ti oc u~^. The figure indicates r„ocu^^asu— >oo 
for all n. FigureHI^b) displays Tn/Ti, and we find the re- 
lation Tn — nTi as u — )■ 00. We observe that T„ oc nu~'^ 
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We observe from this simple calculation that R'^ oc 
nu^ in the limit of u — > 0. The reason is that a wave 
packet of a fermion collides n times with fermions of the 
other spin, while it moves to the opposite side (note that 
almost all probability amplitude of it reaches the oppo- 
site). Therefore the total number of collisions is n^ so 
the reflectance R'^'^ = (n^i?i)/n = nRi. This result is 
consistent with our simulation results by t-DMRG. 

In contrast, in the limit of u — )• oo, it is shown that 
1 — i?^'^ = T^'^ ex u~^ and T^'^ is independent of n. In this 
limit, the particles are almost perfectly reflected by each 
other, so passing to the opposite side occurs only on the 
border between the two clusters. Since the collision on 
the surface occurs just n times, the transmittance T^'^ = 
(nTi)/n — Ti. This result conflicts with our simulation 
with the quantum effects fully considered. 

Therefore the quasi-classical picture, "quasi-classical 
independent collision" , is not valid in the limit of u — ?> oo, 
so a different picture to explain the results for u — >■ cxd is 
needed. Since the particles are fermions, the maximum 
number of particles which share a site is just two. These 
results imply that the number of possible one-to-one col- 
lisions on the surface, which is n^, determines both the 
reflectance in the limit of u — )■ and the transmittance 
in the limit of w — > oo because of dynamically emerging 
quantum many-body effects. 



CONCLUSION 



FIG. 4. (a) Transmittance r„ for u = UAx between 144.! 
and 3276.8 [hrjT^^], and (b) ratio of transmittance Tn/Ti. 



in the limit of strong on-site interaction. 

The numerical results indicate that the sign of u has no 
effect within the computational accuracy. We also calcu- 
late for D = 25, and the results i?„ and r„ are changed 
typically by 10%. For smaller values of D, the wave 
function loses its original form faster as time advances 
because of the effect of the discretization. However the 
relations i?„ ~ nRi and T„ — nTi are still observed. 



Using the time-dependent density matrix renormaliza- 
tion group method and the Fermi-Hubbard model, we 
calculate collision dynamics between two kinds of fermion 
clusters in both cases of weak and strong interaction. Re- 
flectance (transmittance) is in proportion to the particle 
number when the interaction is weak (strong). At least in 
the two extremes, the cluster-cluster collision dynamics 
is mainly dependent on the number of possible one-to-one 
collisions. 
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DISCUSSION 
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We assume a simple quasi-classical dynamics model to 
compare it with the calculation results. The assump- 
tion of the quasi-classical model is that the fermions 
in the clusters are independently localized in each po- 
sition, and the collision dynamics is described by a se- 
ries of one-to-one collisions between two wave packets 
of the independent fermions. Under this assumption of 
"quasi-classical independent collision" , we can calculate 
quasi-classical reflectance R"^ only from the one-to-one 
collision reflectance Ri (by deflnition Rf^ — Ri)- 
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